Dynamic roughening and fluctuations of dipolar chains 
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Nonmagnetic particles in a carrier ferrofluid acquire an effective dipolar moment when placed in 
an external magnetic field. This fact leads them to form chains that will roughen due to Brownian 
motion when the magnetic field is decreased. We study this process through experiments, theory 
and simulations, three methods that agree on the scaling behavior over 5 orders of magnitude. 
The RMS width goes initially as t^^'^ , then as t^^'^ before it saturates. We show how these results 
complement existing results on polymer chains, and how the chain dynamics may be described by 
a recent non-Markovian formulation of anomalous diffusion. 
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Magnetic holes Q, Q, |^ are nonmagnetic micrometre 
sized spheres suspended in a ferrofluid, much larger than 
the magnetites in suspension (nm). In an external mag- 
netic field H, these holes acquire an effective dipolar mo- 
ment equal to the opposite of the dipolar moment of the 
displaced ferrofluid. This is true for any non-magnetic 
material that makes up the spheres. When placed be- 
tween two non-magnetic glass plates the spheres aquire 
interactions that may be fine tuned to produce well de- 
fined separation distances and inter particle forces. We 
have recently obtained the analytic form of these inter- 
actions and verified it experimentally Among the 
many intriguing structures that may be predicted from 
this theory are particle chains. 

Collective Brownian motion in particle chains have 
been studied extensively over the last decade. In 
general this motion is governed by a complex interplay 
of different mechanisms such as hydrodynamic particle- 
particle interactions, particle and fiuid inertia, chain- 
chain interactions, Brownian forces and direct viscous 
drag forces. The magnetic hole chains represent a simple 
prototype system that is designed to eliminate all but 
the essential mechanisms needed to produce non-trivial 
collective behavior. The isolation of these mechanisms, 
which are Brownian fiuctuations, viscous drag and an 
anisotropic interparticle attraction, allows a straightfor- 
ward theoretical treatment of the kinetic roughening pro- 
cess. These predictions are confirmed both by our exper- 
iments and simulations over 5 orders of magnitude in the 
dynamic domain. The size of this domain along with the 
theoretical simplicity of the model appears to allow sig- 
nificantly more conclusive statements both on the scaling 
behavior of the chain motion and on the associated pref- 
actors, than in existing works. 

By studying dynamic roughening of initially straight 
chains we make contact with the theory of kinetic inter- 
face growth processes. It is shown that since the long 
range hydrodynamic interactions are eliminated in our 
setup by the presence of confining interfaces at a dis- 
tance shorter than the chain length, the continuum de- 




Figure 1: Typical dynamic roughening of dipolar chains of 52 
magnetic holes of diameter 4/im after a sudden decrease of 
the external magnetic field. Experiments are shown in gray 
and simulations black and white. 



scription of the chain roughening process conforms to the 
Edwards- Wilkinson equation The model developped 
here is also a discrete generaHzation of a Rouse model 
0, first developped to study polymer dynamics. 

In the end we will show how the Markovian N-particle 
description of the entire chain may be be contracted to 
a non-Markovian description of a single particle in the 
chain. This is done by integrating out all interaction de- 
grees of freedom. What is left is a generalized Langevin 
equation with long term memory. It has recently been 
shown how such an equation may be used to predict 
anomalous diffusion exponents 01, and indeed these ex- 
ponents coincide with our independent predictions and 
measurements. 

The ability of particle chains to change the rheological 
properties of their carrier fiuids has given rise to practi- 
cal applications and designs such as dampers, hydraulic 
valves, clutches and brakes 0. In constant fields (elec- 
tric or magnetic, depending on the nature of the dipoles), 
these chains aggregate laterally due to their ther- 
mal fiuctuations. Therefore, the precise quantification 
and understanding of these fiuctuations along an isolated 
dipolar chain is an important component to understand 
the ag greg ation phenomena in constant fields of MR 
or ER [lU fiuids, as well as in systems of magnetic holes 

We note that also by coating the microspheres with 
bio-active materials, such as streptavidin they may 
be used for the direct manipulation of single strands of 
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DNA. For such application the quantitative control of 
Brownian fluctuations is essential. 

In the experiments monodisperse polystyrene spheres 
01 of diameters a — 3.0 or 4.0 /xm were dispersed in 
a kerosene based ferrofluid p^] of susceptibility x — 0-8 
and viscosity 77 = 6 • 10~^Pa s, inside a glass cell of size 
38 mm x 8 mm x 10 /im. A pair of outside coils pro- 
duced magnetic field strengths up to H = 20 Oe. The 
setup was mounted under an optical microscope with an 
attached video camera recording four frames per second. 
Low volume fractions (< 1%) of microspheres were used 
and chains were grown ^3 by applying a constant field 
of about 7J = 18 Oe parallel to the thin ferrofiuid layer 
for about 20 minutes. The ceh was searched for long 
isolated chains of 30-120 spheres. The field was then re- 
duced to a constant value H in the range 2 Oe < < 10 
Oe while the motion of one long chain was recorded and 
analyzed. One pixel of the video image corresponded to 
0.5/im, and the uncertainty in particle position could be 
reduced to 0.2 /xm by utihzing the intensity profile of the 
pixels showing the particle. The experiments illustrated 
in Fig. n] are challenging in part because this accuracy is 
needed to reveal the scaling behavior of the chains. 

In order to obtain an equation of motion for chains 
we define the lateral displacement of particle i, /li, from 
the initial straight line. A particle pair at a separation 
r and angle 9 to the external magnetic field, experi- 
ences a coupling energy U = ^/(t^(1 — 3cos^ 6')/(47rr'^) 
where a — —VxVL [J, and V is the hole's vol- 
ume and X — 3x/(3 -I- 2%) takes into account the de- 
magnetization factor of spherical particles. The ratio 
of the maximum interaction energy over the thermal 
energy ksT characteristic of the random forces due to 
the molecular motion in the solvent is defined as ^3 
A = {fj.fa^)/{2T:a^kT). In Fig.[l]the initial value A = 866 
is reduced to A = 24. Neglecting all magnetic interac- 
tions but the nearest neighbor ones and performing a 
Taylor expansion of the magnetic interaction potential 
around the minimal energy configuration, a straight line 
with spacing a, the lateral component of the magnetic 
force on sphere i is F*^ — a{hi+i — 2hi + hi-i), with 
a = — 6A • kT/a^. Since the Reynolds 

number in this system is very small (typically Re= 10~^), 
the hydrodynamic forces are linear in the particle veloc- 
ity, and = —Khi where k — i-Krja. 

Newton's second law for the i-th sphere is then 

mil, = Ft' + F," + at) (1) 

where the fiuctuating force C,i{t) is due to the molecu- 
lar nature of the fiuid and gives rise to the Brownian 
motion of the particle. At time scales exceeding the vis- 
cous damping time tm = m/ k = a^p/18r] « 10~^ s, the 
inertial term mhi is negligible. Due to the presence of 
confining plates the inertial motion in the fiuid also de- 
cays on this time scale, so that we can neglect any non- 
Markovian corrections to the above equations as well as 



such corrections in the fiuctuating force [13|. We may 
therefore write (Ci(i)Ci(O)) = 2KkBT6{t)5ij, where the 
prefactor refiects the equipartition of particle kinetic en- 
ergy, i.e. ksT = {mhf) [ijj. Combining the above equa- 
tions Eq. (QJ can be written 

a 1 

h,^ -{h,+i + h,^i-2h,) + -Ciit) . (2) 

K K 

For spatial scales above a and times above irn, the above 
reduces to the Edwards-Wilkinson equation Kdh/dt = 
aa?d'^h/d'^x^ + C,i{t) , also known in polymer dynamics 
as the Rouse model [fi]. 

We consider an isolated chain of N particles, and are 
interested in the dynamic roughening of the chain. To 
observe this experimentally A is decreased from a value 
Ao ^ 1, to a finite value still greater than 1 (to ensure 
that the chain does not melt), and the root mean square 
width of the displacements along the chain is recorded. 

It is convenient to describe this by using the dis- 
crete space-Fourier transform along the chain ft.„ = 
for n = 0...N - 1. Equation Q 

then takes the form 

hn = -(^nK +Cn/K , (3) 

with the dispersion relation — 2a{l~cos{2Trn/N)) / k , 
and random terms obeying {Cm(t)C*{0)) — 
2KkBTS{t)Smn/N . Equation JSJ is easily solved to 
give 

hn{t) = /i„(0)e-""* + I dt' e-'^"(*-*')C„(t')/'« • (4) 

Setting hn (0) — and taking the thermodynamic average 
of the square of the above equation leads to 

{hn{t)h*^{t)) = ^ l_e-^ 

for n^O, and(|/io(i)P) = {2kBT/KN)t. It is seen from 
Eq. 101 that each Fourier mode is initially in a free dif- 
fusion regime, (|/in(i)P) ^ t, for t <^ Tn = ll2ujn and 
saturates when t^Tn — \/2ijJn- The minimum and max- 
imum saturation times are respectively T[^/2 — t — K/%a 
and Ti = TV^t/tt^ for the shortest and longest wave- 
length. 

We are interested in the mean square width of the 
chain = E^ih^ - Ejh,/Nr/N = En=iiK\'). 
When t <^ T all Fourier modes are in free diffusion and 
W"^ ~ 2kBTt/ n . This result arises only because of the 
existence of a shortest wavelength a in the system. In the 
continuum limit a ^ 0, r ^ and this regime does not 
exist. Later on the modes associated with progressively 
longer wavelengths reach their saturated states and this 
is refiected in a new scaling behavior of W . By inserting 
Eq. @ in the above expression for W"^ it is straightfor- 
ward to show that it satisfies the Family- Vicsek scaling 
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Figure 2: Scaling of the random mean square width of the 
transverse displacements: theory, experiments and brownian 
dynamics result. The units are Oe, /im and seconds for H, a, 
W and t respectively. 



Figure 3: Power spectrum of the space Fourier transform of 
the fluctuations normalized by time: scaling of the autocor- 
relation function in time and space. The theory is given by 
Eq. ©. 



form W'^ = NF (t/N'^). Moreover, the exact form of F 
may be obtained. Along with the expressions for r, k 
and a this gives the result 
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when t <^ T 
fit when T < t < tN'^I'k'^ 



when t > tN'^ ItP- 



The hydrodynamic coupling between the particles and 
the confining plates was taken into account by renormal- 
izing the drag coefficient as 77/770 — 2/[l — 9a/16d] — 1 = 
1.40 (1.7) for a = 3nm{-ifim) This scaling law 

was checked in 15 experiments where chains of = 36 
to 59 spheres of diameters 3 or 4 /xm were allowed to 
evolve from states where A was reduced to values between 
2.7 and 267. In Fig. [2l the average of the scaled width 
WHyJajN is displayed as a function of the scaled time 
H'^t/N'^ . The average is taken both over time intervals 
of Q.Qlt and over different experiments. 

The relevance of the long-range nature of the dipo- 
lar interactions and of the linearization of the magnetic 
interactions was studied via Brownian dynamics simula- 
tions, where Eq. Q without the mho term was solved, 
this time with the full dipolar form of the magnetic in- 
teractions computed for every particle pair. A repulsive 
potential ~ exp{—\QQr / a) when r < a was used to pre- 
vent any significant overlaps. The average width over 
100 simulations was evaluated for A^ = 36, a = 3^m and 
H = A Oe. Figure 12 demonstrates that both the sim- 
ulations and theoretical results agree with experimental 
measurements in all three scaling regimes, although the 
crossover between the t^/'^ and t^/^ regimes extend over a 
full decade. The W ^ t^/^ regime is visible over roughly 
2.5 decades. The high variance around the average of 
W in the saturated regime comes from the small number 



(3) of experiments that were carried out at this reduced 
time, along with the fact that this regime is sensitive 
to the longest wavelength, for which the variance of the 
amplitude is highest (it scales as l/w„). 

It is also possible to obtain the dynamic scaling from 
the equilibrium behavior. For this purpose consider the 
departure A/i^ from an arbitrary initial configuration, i.e. 
Aft,i(t,to) = hi{t+to)—hi{to). Using Eq. 1^ again we find 
that the averaged space Fourier transform of A/ii(i,to) 
obeys 



(|Afe.4t,to)p) 
At 



-iJnt 



UJr>t 



(6) 



where A = 2kBT/37rr/aN. Above we have used the satu- 
ration level (|/in(to)P) = ksT / uNuJn from Eq. Com- 
paring this expression with Eq. © shows that the lateral 
fiuctuations at equilibrium behave similarly to the lateral 
displacements starting from a straight chain during the 
non-equilibrium roughening stage, but display an ampli- 
tude times larger and develop 2 times more slowly. 

In Fig. El we have compared the theory of Eq. (EJ to 
A^ = 57, a = i^m experiments with A = 113 and 18.1, 
for which r = 0.022 s and 0.14 s. For each t the power 
spectrum |A/i„(t, io)P, calculated using a Hamming win- 
dow, was averaged over all possible t^s in the 30 minutes 
the experiment lasted. 

The agreement between these experiments and this 
theory, where there are no free parameters, is satisfactory. 
Note that w„ was evaluated using its discrete form. Us- 
ing the continuous asymptotic limit of the Rouse model, 
uon — LOiin? instead, would result in systematic deviations 
from this model for large wavenumbers. 

Finally, we sketch the structure of the connection be- 
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tween the present problem and the general framework of 
anomalous diffusion of particles in systems with memory 
We consider particle z = and use Eq. Q to inte- 
grate out the i ^ variables. By applying time Fourier 
transforms and neglecting all but the z = mass it is 
possible to arrive at 

mho{uj) = -mr{uj)ho{uj) + F{uj) (7) 

where the time Fourier transform is denoted by hat 
symbols, F(t) is the fluctuating part of the force and 
r(cj) = [A / m) aK / ioj for 1/ri <^ lo <^ 1/t and ^{lu) ^ 
{A/m)yJaKTi/i for lu <^ 1/ti. The inverse Fourier trans- 
form of Eq. is a Generalized Langevin Equation of 
the Mori-Lee form 01, 

■mho{t) = -mf dtiT{t-ti)ha{ti)+F{t). (8) 

From the exact form of F{t) the fluctuation dissipa- 
tion theorem for such a generahzed Langevin equation, 
{F{t)F{0)) = 2m/csTr(|i|), can be checked directly. The 
memory effects and correlations in F{t) come from the 
interaction between the rest of the chain and the fluid. 
According to theory 0, the random mean square width 
of the observed particle's displacement scales at long time 
as W"^ ~ f-" if f(w) ~ in the limit w ~ 0. 

In the present case, the asymptotic behavior of the re- 
sponse function corresponds to a = for flnite chains 
and t > Ti, and to a = 1/2 for times below ti (ti oo 
when N —> oo). The observed behavior ^ t^^^ for 
t <C Ti is then in agreement with this theory, and the 
present system is a simple experimental example validat- 
ing this theory. 

For chains in an unbounded medium, the Zimm model 
predicts W ~ t^^^ @|, and experiments carried on vari- 
ous MR fluids report W ~ ^o.asio.os ^ r^j^^g behavior 
is attributed to either hydrodynamic- or long rang dipo- 
lar interactions. In our system hydroynamic interactions 
are suppressed by the walls. This may be understood 
by noting that the flow in our system is described by 
a simple Darcy law at length scales above d. This im- 
plies that a local flow perturbation decays as 1 /r^ which 
is faster than the 1/r decayjgiven by the Oseen tensor 
in an unbounded medium Since the chain length 
is typically much larger than d, hydrodynamic interac- 
tions are local and do not affect the longer wavelengths 
of our chain motion. On the other hand we note that if 
our experimental and simulation data were truncated as 
logiQ{H^t/N^) < -1 they could support the W ~ t^/s 
interpretation made by Furst & Gast fl. As these au- 
thors attribute their exponent 3/8 to long range dipolar 
interactions, our results may indicate that such interac- 
tions do play a role, but only in a crossover regime to the 
W ~ i^/^ behavior. 

In conclusion, we have established the scaling behav- 
ior of conflned dipolar chains over 5 orders of magni- 



tude by matching experiments, theory and Brownian dy- 
namics simulations. In passing we have made contact 
with nearby theories of both polymer models and non- 
Markovian formulations of anomalous diffusion. These 
results are attractive in particular because of the versa- 
tility and easy control of the experiments. 
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